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The paper orders certain important issues related to both uncorrelated and correlated networks 
with hidden variables. In particular, we show that networks being uncorrelated at the hidden 
level are also lacking in correlations between node degrees. The observation supported by the 
depoissonization idea allows to extract distribution of hidden variables from a given node degree 
distribution. It completes the algorithm for generating uncorrelated networks that was suggested 
by other authors. In this paper we also carefully analyze the interplay between hidden attributes 
and node degrees. We show how to extract hidden correlations from degree correlations. Our 
derivations provide mathematical background for the algorithm for generating correlated networks 
that was proposed by Boguna and Pastor-Satorras. 
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I. INTRODUCTION 

Recently, the techniques of equilibrium and non- 
equilibrium statistical physics were developed to study 
complex networks 0, 0, @- This paper is devoted 
to equilibrium correlated networks, with a special at- 
tention given to networks with two point correlations 

inula. 

What does it mean that a network is correlated? In 
simple words one can say that there exist certain rela- 
tionships between network nodes. For example, when 
one considers a social network i.e. a group of people 
with links given by acquaintance ties one may expect 
that young people are mostly surrounded by other young 
people. One may also expect that wealthy individuals 
are more often associated with other wealthy individuals 
than with poor once. In some sense, the above exam- 
ples let one suppose that social networks are positively 
correlated, at least when one considers individual's age 
or income. The situation is more contentious if one asks 
for relationship between gender of the nearest neighbors. 
Now, it is difficult to assess if social networks are posi- 
tively or negatively correlated. The above examples show 
that even taking into account a single network one can 
observe correlations at different levels of the system com- 
plexity. Each node i in such a system has assigned a set 
of different attributes like: gender gi, age a^, education 
e^, attractiveness ki etc. The last property may be quan- 
tified as a number of nearest neighbors of the considered 
individual. In the graph theory |9j the quantity is known 
as node's degree. 

The above outlined network correlations and multilevel 
structure constitute two main issues undertaken in this 
paper. 

Multilevel topology has been recently considered by 
several groups of researchers [l(J, HH [I2> LLJ 03 an d, at 
present, the proposed models are known as random net- 
works with hidden variables. To be precise, to our knowl- 
edge, none of the proposed models considers the number 



of complexity levels larger than two |jj i.e. each node 
i is characterized by at most two parameters: hidden 
attribute hi and node's degree ki. In general, random 
networks with hidden variables have a fixed number of 
vertices N . Each node in a network belonging to this 
class of models has assigned a hidden variable hi (fitness, 
tag) randomly drawn from a fixed probability distribu- 
tion R(h) (throughout the paper we use the symbol R 
with reference to distributions at hidden level and P with 
reference to node degrees) . Edges are assigned to pairs of 
vertices {«, j} with a given connection probability ry. In 
the simplest case ry depends only on values of the hid- 
den variables hi and hj , but in a more general situation it 
may for example involve hidden variables characterizing 
the nearest or the next-nearest neighbors of both con- 
sidered nodes i and j. In fact, the first case represents 
networks with at most two-point correlations at the level 
of hidden variables, while the latter one allows for higher 
order correlations. 

Below we outline the concept of network correlations 
in a more rigorous way. The introduced ideas will be 
completed and widely discussed in the next section. 

From the mathematical point of view, the lack of two 
point correlations means that the probability R(hj/hi), 
that an edge departing from a vertex i of property hi ar- 
rives at a vert ex j of property hj , is independent of the 
initial vertex i |lall6| . The above translates into the fact 
that the nearest neighborhood of each node is the same 
(in statistical terms). On the other hand, when R(hj/hi) 
depends on both hi and hj one says that the studied net- 
work has two-point correlations 0, • To characterize 
this type of correlations one usually takes advantage of 
the joint, two-dimensional distribution R(hi, hj) that de- 
scribes the probability of randomly chosen edge to con- 
nect vertices labeled as hi and hj. 

In order to characterize network in a more de- 
tailed way the concept of higher order correlations, 
given by multidimensional probability distributions, 
should be used. In this paper we limit ourselves 
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to two-point correlations. The lack of higher or- 
der correlations is ensured by the factorization of 
the conditional distribution fl(/iW,ft( 2 ',...,li"/li) = 
R(hW/h)R(hW/h)...R(hW/h), that describes the 
probability of a node of property h to have x neighbors 
labeled as h^ 2 \ . . . , h^ x \ Such networks, with the 
only two-point correlations are called Markovian, due to 
the reason that they are completely defined by the joint 
distribution R(hi, hj), which on its own turn completely 
determines the conditional distribution R(hj/hi) and the 
distribution of hidden variables R(h). Relationships be- 
tween R(hi, hj), R(hj/hi), Tij and R(h) will be analyzed 
later. 

After this short introduction to the concept of hidden 
variables and to the problem of network correlations we 
may go back to the main topic of the paper i.e. the 
interplay between different levels of the network com- 
lcxity. The approach has been initiated by Soderberg 
and developed by Boguha and Pastor-Satorras 0. 
Boguha and Pastor-Satorras have concentrated on the 
question: how correlations at the level of hidden vari- 
ables affect pattern of connections at the level of node 
degrees. Given R(hi, hj) the authors have derived an an- 
alytical expression for the joint distribution of degrees of 
the nearest neighbors P(ki,kj). In this paper, we ask 
the reverse question: what kind of hidden correlations 
R{hi, hj) produces the given pattern of node degree cor- 
relations P(ki, kj). 

As a matter of fact, since most of us are better ac- 
quainted with node degree notation than with abstract 
hidden variables, the reverse approach seems to be very 
interesting, at least from the methodological point of 
view. It is already well-known that there exist degree cor- 
relations in real networks 4] . On the other hand, due to 
the lack of data nothing is known about correlation at the 
hidden level, from which the observed network structure 
arises. The paper represents a small step towards un- 
derstanding the phenomenon of self-organization in com- 
plex networks beyond the predominant approach of the 
so-called evolving networks 0, . 

The paper is organized as follows. In Sec. [n] we re- 
view general results on correlated random networks with 
hidden variables 7]. The salient issues concerning two- 
point correlations are discussed in this section. Sec. IIIII 
is devoted to theoretical aspects of our reverse problem. 
Some remarks on uncorrelated networks and a practical 
algorithm of generating random networks with two-point 
degree correlations are also given in this section. Finally, 
in Sec. IIVI we draw our conclusions. 



II. 



CORRELATED NETWORKS WITH HIDDEN 
VARIABLES 

A. Construction procedure 



Probably the most attractive feature of random net- 
works with hidden variables is the construction procedure 



that consists of only two steps: 

i. first, prepare N nodes and assign them hidden 
variables independently drawn from the probability 
distribution R{h), 

ii. next, each pair of nodes link with a probabil- 
ity nj . 

Note that the above procedure does not include any am- 
biguous instructions like in the case of . . . do .... The 
lack of ambiguities enables analytical treatment of the 
discussed models. 

The particularly simple model of sparse random net- 
works with hidden variables has been recently considered 
by Boguha and Pastor-Satorras. Provided that the con- 
nection probability scales according to 



f(hj,hj) 
N 



(1) 



where j(hi, hj) = f(hj, hi), the authors have shown that 
the degree distribution of nodes characterized by a fixed 
value of hidden variable h is given by the Poisson distri- 
bution 



g(k/h) 



-huh 



k\ 



(2) 



The last result is very important because it joins the 
two levels of the network complexity and allows to freely 
move between the notation of hidden variables and the 
notation of node degrees (see section ITTT1 . In particular, 
taking advantage of one can see that the average 
connectivity of nodes characterized by h is simply h 



(k){h) = Y. k 9( k l h ) = h ' 



(3) 



i.e. hidden variables play the role of desirable degrees. 
Next, due to J2Jl, the degree distribution P{k) gains a 
simple form of convolution 



P(k) = / R(h)g{k/h)dh. 



(4) 



The last relation between both distributions P{k) and 
R(h) implies relation between their moments 



(h n ) = (fc(k-l)...(fc-n+l)>, 



and respectively 

(h) = <fc), 



(h 2 ) = (k(k-l)). 



(5) 



(6) 



As a matter of fact, in the case of other pair-related dis- 
tributions, like R{hi,hj) and P(ki,kj), it is also possible 
to obtain relations similar to (QJ. However, before we 
proceed with the interplay between hidden variables and 
node degrees, we would like to discuss the issue of two- 
point correlations. 
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FIG. 1: Schematic representation of an ensemble of random 
networks with a given two-point correlations at the level of 
hidden variables. Detailed description of the figure is given in 
the text. 



B. Two-point correlations 

One usually thinks about a network as about a col- 
lection of nodes and links. Distribution of hidden vari- 
ables R(h) is the basic characteristic of the set of nodes, 
whereas the joint distribution R(hi,hj), describing two- 
point correlations, applies to the set of edges (see FigQ. 
Both distributions express the probability that a random 
representative of its own ensemble has assigned a given 
attribute i.e. property h in the case of node and pair of 
hidden variables {hi, hj} in the case of edge. 

At the moment, let us briefly comment on the Fig^ 
The second stage of the figure shows the set of nodes and 
the set of links, both corresponding to the simple network 
presented at the first stage of the same figure. Obviously, 
basing only on the two sets is almost impossible to recre- 
ate the original network. Such a representation related 
to a single network neglects much information. What is 
lacking are higher order correlations. On the other hand, 
the joint distribution R{hi,hj) characterizing the set of 
edges conveys all the information that is required to con- 
struct an ensemble of Markovian random networks (see 
the third stage of the Fig^]). In fact, all the calculations 
performed in this paper are related to such an ensemble. 

For a given network in the thermodynamic limit, both 
distributions R(h) and R(hi,hj) may be defined in the 
following way 



R(hi, hj) 



E{hj, hj) 
E ' 



(8) 



where N(h) gives the number of nodes assigned by h and 
N = J N(h)dh, whereas E(hi,hj) describes the number 
of links labeled as {hi, hj} and due to the normalization 
condition E = J J E(hi, hj)dhidhj. 

The connection probability r»j Q may be simply ex- 
pressed by the ratio between the actual number of edges 
connecting vertices of tags hi and hj and the maximum 
value of this quantity 



E(hj, hj) 

Ejnax{hi, hj) 



(9) 



Since, during the network construction process one an- 
alyzes all pairs of nodes, the maximum number of the 
considered connections is given by E max — N{hj)N{hj). 
Now, taking advantage of (JHJ) the expression for r%j may 
by rewritten in the following form 



{h}R(hj,hj) 
Flj NR(hi)R(hjY 



(10) 



To fulfill discussion on two point correlations it is nec- 
essary to find relation between the two distributions (JJJ 
and (JBJ ■ The relation is encoded in the so-called detailed 
balance condition 



R(hj/hi) 



(fe)fl(/ti,/tj) 
hiR(hi) 



(11) 



Summing both sides of the last identity over hj 
(J R(hj/hi)dhj = 1) one obtains 



R(hi 



(/,) 



R(hi, hj)dhj 



(12) 



As a truth, almost all the above derivations that have 
been done for hidden variables also hold for node degrees 
(see [Tfll| ^1 . In particular, the degree detailed balance con- 
dition possesses the same form as (ITTTl 



(k)P(ki,kj) 



(13) 



and respectively the degree distribution in Markovian 
networks is given by the relation analogous to l|12l) 



P(ki)= { -^J2 P ^^ 



(14) 



Interplay between hidden variables and node 
degrees 



(7) 



Each node of the considered networks is characterized 
by two parameters: hidden variable h and degree k. The 
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probability p(k n h) that two given parameters h and k 
meet together in a certain node is described by the below 
identity 

p(k Dh) = R{h)g(k/h) = P(k)g*(h/k). (15) 

The meaning of both conditional distributions g{k/h) 
and g*(h/k) is simple. The first distribution g(k/h) 
has just been introduced. It describes the probability 
that a node labeled by h has k nearest neighbors. The 
second distribution g*(h/k) is complementary to the first 
one and gives the probability that a node with k nearest 
neighbors is labeled by h. 

The knowledge of g(k/h) and g*(h/k) allows one to 
find the relation between the joint distributions R(hi, hj) 
and P{ki, kj) characterizing pair correlations respectively 
at the level of hidden variables and at the level of node 
degrees. Simple arguments let one describe the interplay 
by the following relation po| 

P(ki,kj) = J J g(ki — l/hi)R(hi,hj)g(kj — l/hj)dhidhj. 

. (16) 

The last expression states that if one knows hidden 
correlations then it is possible to calculate degree cor- 
relations. Relating the problem to the present state of 
the art in the field of complex networks it is a very ar- 
tificial situation. The concept of networks with hidden 
variables is still little-known because most of researchers 
are strongly attached to the node degree notation. From 
this point of view, the reverse problem that would answer 
the question: what kind of hidden correlations produces 
given degree correlations, seems to be very attractive. 
In addition, since one knows how to construct Marko- 
vian networks with hidden variables, the solution of the 
above problem would also allow, in a simple way, to gen- 
erate networks with given two-point degree correlations. 
We work out the problem in the next section entitled 
" Degree correlations from hidden correlations" . 

In fact, the idea of generating networks with given de- 
gree correlations by means of networks with hidden vari- 
ables originates from Boguha and Pastor-Satorras. The 
authors have argued that since the conditional proba- 
bility g(h/k) is Poissonian (J2J, the joint distribution for 
node degrees Hlfi|) approaches the joint distribution for 
hidden variables for ki , > 1 



P{ki, kj) — R{ki, kj), 



(17) 



and respectively asymptotic behavior of the degree dis- 
tribution is given by 



III. HIDDEN CORRELATIONS FROM DEGREE 
CORRELATIONS 

A. Depoissonization 

At the moment, let us examine Eq. Q more carefully. 
Depending on whether h is discrete or continuous vari- 
able the expression for P(k) is respectively discrete or 
integral transform with Poissonian kernel |21| . The issue 
of determining R(h) from P(k) is simply the problem of 
finding the inverse transform 22] and one can show that 
for a given P(k) there exists a unique R(h) which satisfies 
When h is continuous one gets (see Appendix) 



R(h) = e h T- l [G{ix)] 



(19) 



where T 1 denotes the inverse Fourier transform and 
G(x) describes generating function for the degree distri- 
bution P(k) 



G{x) = $> fe P(fe). 



(20) 



In the case of discrete h, the inverse Poisson transform is 
given by the formula 



R{h) = e' l Z- l [G(-\ux)}, 



(21) 



where Z 1 describes the inverse Z-transform. 

The same applies for the joint degree distributions 
characterizing two-point correlations. Since 



g(k - 1/h) = ~ g(k/h), 
the formula l|16|l may be rewritten as 



(22) 



''(*■'■■ /■•.;> _ f f g{kj/hj) hj) g{k . /h . )dh . dh ,_ 



hihj 



(23) 

Now, it is easy to see that P(ki,kj) is connected with 
R(hi, hj) by the two-dimensional Poisson transform and 

= e h i+hj jr-l {J r-l [G * {iX)m (24) 



hihj 



where 



P{ki) ~ R{h). 



(18) 



P{k>ii kj 
ki kj 



(25) 



The range of convergence of the two distributions has 
been estimated by Dorogovtsev who has shown that 
both approximations (|17|) and i|18|) are only acceptable 
when R{hi, hj) and R(h) are sufficiently slowly decreas- 
ing. 



The last two expressions can be simply translated for the 
case of discrete hidden variables. Since however Fourier 
transforms are more convenient to work with then Z- 
transforms, we have decided to pass over such a reformu- 
lation. 
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Some remarks on uncorrelated networks 



where 



As we said at the beginning of this paper - the lack of 
pair correlations at the level of hidden variables means 
that the conditional probability R {hj/hi) does not de- 
pend on hi (to differentiate between uncorrelated and 
correlated networks the characteristics related to the for- 
mer case have been denoted by the subscript V). In fact, 
it is simple to show that 



Ro(h 3 /h t ) = ji-Roihj), 



(26) 



and respectively the joint distribution gains a fac- 
torized form 



R {hi, hj 



hjR (hi)hjR (hj) 

W ' 



(27) 



Inserting the above expression into © one gets the for- 
mula for the connection probability in uncorrelated net- 
works 



r °.(h h-)= hihj 



(28) 



Now, the question is: does the lack of pair correlations 
at the hidden level translates to the lack of pair correla- 
tions between degrees of the nearest neighbors? Inserting 
(|27|l into 1|16|) . and then taking advantage of the degree 
detailed balance condition (|13[) one gets the answer 



(k) 



Po(kj), 



(29) 



i.e. the lack of hidden correlations results in the lack of 
node degree correlations. Thus, in order to generate un- 
correlated network with a given degree distribution P(k) 
one has to: 

i. prepare the desired number of nodes N, 

ii. label each node by a hidden variable randomly 
taken from the distribution R(h) given by l(T§|) or 

iii. each pair of nodes link with the probability r°j (|28p. 

(See also other methods for generating random uncorre- 
lated networks with a given degree distribution 

EE3). 

Note also that, due to properties of the Poisson prop- 
agator @ 



(k{hi)) and hj = (k(hj)), 



(30) 



therefore the connection probability r°j(hi,hj) l|28|) may 
be considered as the connection probability p°j(ki, kj) be- 
tween two nodes of degrees ki and kj averaged over all 
pairs of nodes possessing hidden variables respectively 
equal to hi and hj 



o /j i \ kikj 



(32) 



The last expression describing connection probability in 
uncorrelated sparse networks has been already used by 
several authors (in particular see [25ll26l I27l l2q ) . 



C. Examples of uncorrelated networks 

1. Classical random graphs of Erdds and Renyi. The 
degree distribution in ER model is Poissonian 



P(k) 



-<fc> 



fc! 



k > 0. 



(33) 



The first step towards calculation of the required distri- 
bution R(h) is finding characteristic function for Poisson 
distribution. Taking advantage of (|20|l one gets 



G(x) 



,{k){x-l) 



(34) 



Now, inserting the exponential function into (|19fl or 121|) 
one can see that the distribution of hidden variables 
in classical random graphs is respectively given by the 
Dirac's delta function (in the case of continuous h) or 
the Kronecker delta (in the case of discrete h) 



R{h) = 8(h-(k)). 



(35) 



The above result means that all vertices are equivalent 
and the connection probability at the level of hidden vari- 
ables (P%)l is given by 



(fc) 
N ' 



(36) 



Fig. |3 presents probability p°j(ki,kj) of a connection 
between two nodes characterized by degrees ki and kj 
in ER model and other uncorrelated networks. As one 
can see, there exists a very good agreement between the 
formula (l.'S2l) and results of numerical calculation. 

2. Networks with P(k) given by exponential distribu- 
tion. Now, let us suppose that 



P(fc) = 



{k) k 



\ + {k)Y 



fe+i 



k > o. 



(37) 



The generating function of the above degree distribution 
is given by the sum of the geometric series 



G(x) 



1 



1 + (k) - (k}x 



(38) 



and l|19|) provides the exponential distribution of hidden 
variables (see Fig. |3J) 



R{h) 



,-h/(k) 



h>0. 



(39) 
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FIG. 2: Connection probability p°j(ki, kj) in classical random 
graphs (A), uncorrelated networks with exponential degree 
distribution 1371 (B) and asymptotic scale-free networks (1441) 
(C). Points represent results of numerical calculations aver- 
aged over 500 networks of size N = 5000. 



3. Scale-free networks. In mathematical terms, the 
scale- free property translates into a power-law degree dis- 
tribution 



P{k) = —, k>l (40) 
where a is a characteristic exponent and A represents 
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FIG. 3: Distribution of hidden variables and degree distri- 
butions obtained from numerical simulations of uncorrelated 
and correlated exponential networks. 



normalization constant. Generating function for this dis- 
tribution is given by the polylogarithm 



k—l 



(41) 



To derive distribution of hidden variables that leads to 
uncorrelated scale-free network one has to find the inverse 
Fourier transform of the polylogarithm with imaginary 
argument 



R(h) = e h J r - 1 [Li 7 (ix)], 



(42) 



or the adequate inverse Z-transform. 

Unfortunately, it does not appear that a closed-form 
solutions for both inverse transforms can be simply ob- 
tained. Nevertheless, some asymptotic results for scale- 
free networks can be derived. In particular, one can show 
that power law distribution of hidden variables 



R(h) 



(a - l)m( a - 1 '> 



h > 



(43) 



leads to asymptotic scale-free networks with degree dis- 
tribution given by 



P(k) = (a-l)m (a - 1] 



T{k 



a 



l,m) 



fc! 



(44) 



where T(x, y) is incomplete gamma function. In the limit 
of large degrees k 3> 1 the above degree distribution de- 
cays as P(k) ~ k~ a (see Fig. QJ. The effect of structural 
cutoffs in power law distributions of hidden variables with 
a < 3 l)43p. imposing the largest hidden variable to scale 
as hmnx, ~ VN (the relation follows from r,° < 1 l)55)l) 
[T^. E4I l2flj. does not represent any problem in the stud- 
ied formalism. The effect of h max in the scale-free R(h) 
may be considered as an exponential cut-off 



R(h) 



(a - l)m^- 1 '> 



■ exp 



(45) 
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ii. next, label each node by a hidden variable ran- 
domly taken from the distribution R(h) 1191) . where 
P(k) is is given by l|14fl. 

iii. finally, link each pair of nodes with the probability 
rij @, where R(hi,hj) is calculated according to 
the formula lf23|) . 

Although very clear, the above procedure suffers a 
certain inconvenience: given a joint degree distribution 
P(ki, kj) the closed form solution for R(hi, hj) (|24|l is of- 
ten hard to get. Since however there exists a number of 
algorithms for numerical inversion of Fourier transform 
the above does not represent a real problem. 



FIG. 4: Degree distributions in uncorrelated networks with 
scale- free distributions of hidden variables (1431 . Points repre- 
sent results of numerical calculations whereas solid lines cor- 
respond to Eq. lEHl . 



Due to properties of the Poisson transform [21| , the above 
R(h) results in truncated degree distribution (|4*4"|) 

m = T(k-a + l,m) («-l)m^ 



fc! 



As expected (see discussion of Eq. $l<i\i ). in the limit of 
large degrees the last formula approaches iTEST) . 

Before we finish with scale-free networks let us, once 
again, concentrate on the Fig. The figure presents 
degree distributions P(k) in sparse networks with scale- 
free distributions of hidden attributes R(h) |@3J. One 
can observe that although both distributions converge in 
the limit of large degrees, there exist serious deviations 
between (|43H and 144|) in the limit of small degrees. The 
relative behavior of the two distributions let one expect 
that the correct R(h) reproducing the pure scalc-frec de- 
gree distribution P(k) should describe a kind of conden- 
sate with a huge number of nodes characterized by a very 
low values of hidden attributes. On the other hand, de- 
spite the ambiguous behavior of P{k) for small degrees, 
the power-law tail is interesting on his own right. A num- 
ber of real networks have fat tailed degree distributions. 
The above allows us to deduce on fat tailed distributions 
of underlying hidden attributes assigned to individuals 
co-creating the considered systems. 



D. How to generate correlated networks with a 
given degree correlations 

Here, we again make use of the depoissonization idea 
proposed in Sec. IIII Al The procedure of generating ran- 
dom networks with a given two-point, degree correlations 
P(ki,kj) is as follows: 

i. first, prepare N nodes, 



E. Examples of correlated networks 

To make derivations of previous sections more con- 
crete, we should immediately introduce some examples 
of correlated networks. In order to simplify the task we 
will take advantage of general patterns for joint degree 
distributions with two-point assortative (a) and disassor- 
tative (d) correlations that was proposed by Newman [f| 

P d {k ll kj) = P(ki)Q(kj) + Q(h)P(kj) - Q(ki)Q{kj) (47) 

and 

P a (k tl kj) = 2P(ki)P(kj) - P d (h, kj), (48) 

where P(k) and Q[k) are arbitrary distributions such 

that J2kP(k) — T^fe Q(k) = 1- Assortativity / corre- 
lation coefficients [4J for the above distributions are re- 
spectively equal to 



and 



r = — r 



(49) 



where /i p represents the expectation value for P(k), \i q 
has the same meaning for Qlk), whereas corresponds 

to the variance of P(k). Now, in order to facilitate further 
calculations let us assume that 

^( fc ) = 77X P ( fc ) and to = A^( fc )' ( 5 °) 

(q) 



where 



(k)=J2 kp ( k ) and {a)=^kQ{k). (51) 



Note that there is no conflict of notation in the last as- 
signment. Putting the two expressions l|50Jl into P a or 
P d and then taking advantage of the degree detailed bal- 
ance condition (|14f) one can easily check that P(k) corre- 
sponds to degree distribution. Now, given P{k) one can 
execute the first two steps (i. and ii.) of the construction 
procedure described in the previous section. 
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Inserting the relations l|50|) into (|47|l and (|48fl one ob- 
tains (compare it with (|23[l ) 



P d {h,kj) _P(h)Q{kj) Q(h)P(kj) Q(ki)Q(kj 



(k)(q) 



<*><?} 



(?> 5 



52) 



and 



P a {k u k 3 ) =2 P(k i )P(k j ) P d (h,kj) 



(53) 



Due to linearity of the Poisson transform the joint distri- 
butions l|52|l and l|53|l turn out to be particularly useful 
for our purposes. The usefulness means that whenever 
closed form solutions for the inverse Poisson transforms 
of P(k) and Q{k) exist, one can also obtain closed form 
solutions for the joint hidden distributions R d (hi, hj) and 

R a {h u hj) 121. 

Generating functions G*{x,y) for and 
are respectively given by 

G * (x y ) = G p( x ) G i(y) + G q {x)G p {y) _ G q {x)G q {y^ 



(k)(q) 



(*)(?) 



and 



where 



Gt(x,y) = 2 G ^/ y) -G* d ( X> y), (55) 



G p (x) = ^VP(fc) and G q (x) = ^VQ(fc). (56) 

ft k 

Although visually quite complicated all the above formu- 
las are in fact very simple. Now, in order to perform the 
last step (iii.) of our procedure aiming at constructing 
correlated networks, one has to calculate the joint distri- 
bution R(hi, hj). Taking advantage of (|24Jl one gets 

i? d (/t i7 hj) = R(hi)S{h 3 ) + S{hi)R{hj) - S{hi)S{hj) 

(57) 



and 



R a (hi, hj) = 2R(hi)R(hj) - R d {h t , hj), (58) 
where, similarly to Ij50|l one has 



h 

w 



R(h) = -^jR(h) and R(h) = e h T- 1 [G p (ix)], (59) 



and also 



h 

w 



S{h) = -^S(h) and S{h)=e h T- 1 [G q {ix)]. (60) 



Note, that R(h) given in (|59|l expresses distribution of 
hidden variables in the considered correlated networks 
(i.e. the inverse Poisson transform of the degree distri- 
bution P(k)). 



V 



Exponential networks 

■ assortative networks r—0.11 

A uncorrelated networks r-0 

• disassortative networks r--0.ll 




FIG. 5: Average nearest neighbor degree {k nn )(k) for assor- 
tative, uncorrelated and disassortative networks with expo- 
nential degree distributions (see section lill 



Now, let us translate the general considerations into 
a specific example. Suppose that we are interested in 
networks with exponential degree distribution (|37ll 



P(k) 



(k) k 



(1 + <fc)) fe +! 



k > 0. 



(61) 



Distribution of hidden variables in such networks is given 
by 



R(h) = 



-h/(k) 



h>0. 



(62) 



For mathematical simplicity, let us assume that the dis- 
tribution Q(k) responsible for correlations is also expo- 
nential 



Q(k) 



(1 + <<7» fe +i 



k > 0. 



(63) 



Given P(k) and Q(k) one has to ensure that the joint 
hidden distributions l|57l) and l|58l) and also the connec- 
tion probability are positive and smaller than 1. It 
is easy to check, that in our case the condition translates 
into the relation 



V2-(k)- 



(64) 



Now, let us briefly examine the role of Q(k). Depend- 
ing on the value of (q) one obtains stronger or weaker 
correlations characterized by l|49l) 



d 2((fc)-(<z» 2 



(k)((k)+l) 



and 



(65) 



In particular, given (k) the expression <|64[l provides pos- 
sible values for the correlation coefficients that may be 
reproduced in the considered networks 



< r < 



(k) 



(k) + l 



(V2-1) 



(66) 
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where r = \r d \ = \r a \. 

In order to check the validity of our derivations we 
have performed numerical simulations of correlated net- 
works (I47|l and i|48|) with partial distributions given by 
(|61|1 and (|63|l . Simulations were done for networks ( both 
assortative and disassortative) of size N = 10 4 , (k) = 2 
and (q) = \[2. Given (fc) = 2, the value of {q) = \/2 
corresponds to the maximum value of the correlation co- 
efficient |r maa; | =0.11 ijnni)- At Fig. |3]we depict results 
corresponding to degree distributions (|61[) in the consid- 
ered networks. Fig. [3] presents effects of node degree 
correlations expressed by the average degree of the near- 
est neighbor (ANND) that is defined as 



{k nn )(k) = Y,hP(h/k). 



(67) 



our findings complete the algorithm for generating ran- 
dom uncorrelated networks that was suggested by other 
authors 0, |53 ■ We also show that the connection prob- 
ability in sparse uncorrelated networks is factorized func- 
tion of node degrees l{3"21) . 

In this paper we also carefully analyze the interplay 
between hidden attributes and node degrees. We show 
how to extract hidden correlations from degree correla- 
tions and how to freely move between the two levels of 
the networks complexity. Our derivations provide math- 
ematical background for the algorithm for generating 
correlated networks that was proposed by Boguna and 
Pastor-Satorras 0. 
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It is already well-known that in the case of uncorrelated 
networks l|29l) the ANND does not depend on k 



<O(*0 = 



(k 2 



(k) 



(68) 



whereas in the case of assortatively (disassortatively) cor- 
related systems it is an increasing (decreasing) function 
of k. One can find that in our example the corresponding 
functions are given by the below formulas 



(Kn)(k) 

(kL)(k) 

(Kn)(k) 



l + 2(fc), 
(k° n )(k)-f(k), 
(k° nn )(k) + f(k), 



(69) 
(70) 
(71) 



where 



f(k) = 2((k)-(q)) 



M 

(k) 



fc-1 



(k) 



k+l 



(72) 

As one can see (Fig. |5J), the fit between computer simu- 
lations and the above analytical expressions is very good, 
certifying the validity of the proposed algorithm (section 
for generating random networks with a given de- 



imp 

gree correlations 



IV. CONCLUSIONS 

In this paper we refer to the set of articles devoted 
to the so-called random networks with hidden variables. 
Importance of this paper consists in ordering certain sig- 
nificant issues related to both uncorrelated and corre- 
lated networks. In particular, we show that networks 
being uncorrelated at the hidden level are also lacking in 
correlations between node degrees. The observation sup- 
ported by the depoissonization idea fsection lill All allows 
to extract distribution of hidden variables from a given 
node degree distribution. Until now the distribution of 
hidden variables required for generation of a given de- 
gree sequence had to be guessed. From this point of view 
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VI. APPENDIX 

The inverse Poisson transform for the case of continu- 
ous h (tT§j) h as been derived by E. Wolf and C.L. Mehta 
in 1964 |22|. Below we outline the derivation, and then 
adopt it for the case of discrete h lj?T|) . 

Our aim it to inverse the formula 



P(k) 



-R(h)dh. 



Let 



F(x) 



i ixh R(h)e~ h dh 



(73) 



E 

fe=0 



(ixh) k 



k\ 



R(h)e-' l dh 



(74) 



X>) fc P(*0 = G(i 



k=0 



where G <|20|) represents generating function for the de- 
gree distribution P(k). Then by the Fourier inversion 
formula one gets the expression i|19|) 



R(h) =e-' l T- 1 [G(ix 



(75) 



The above derivation can be simply adopted for the 
case of discrete transform with Poissonian kernel i.e. 



P(k) = £ 



- ~- h h k 



k\ 



R(h). 



(76) 
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Let us introduce an auxiliary sequence 

J(h) =e-' l R(h). (77) 

It is easy to show that the Z-transform of this sequence 
equals to generating function of the degree distribution 

*\m = £ J -f ( 78 ) 

h=0 

_£(f : t**j}!\ m (79) 

h=0 \l=0 ' ) 



= ^2(lns- 1 ) l P(l) = G(lns- 1 ). (80) 

1=0 

Applying the inverse Z-transform to the last expression 
one obtains the formula 121|) describing distribution of 
hidden variables 

R(h) = e h Z- 1 [G{\ns- 1 )]. (81) 
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